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I .  Introduction 


This  report  describes  a  two  dimensional  stochastic  process 


and  a  program  for  its  simulation.  The  process  can  be  considered 
to  describe  a  random  track  of  predetermined  length  that  origi¬ 
nates  and  terminates  at  predetermined  end  points.  The  process 
is  a  modification  of  a  process  that  is  without  the  length  and 
termination  constraints. >,  This  process,  which  is  also  described, 
is  a  generalization  of  the  random  tour  process  described  by 
Washburn  in  Reference  1. 

A  particular  process  is  defined  by  the  track  length,  the 
maximum  expected  track  segment  length  and  the  distance  between 
the  end  points.  For  motion  at  a  constant  speed,  the  track 


length  constraint  is  equivalent  to  a  transit  time  constraint. 

a 

In  Section  IT,  the  process  and  the  unmodified  process  are 


described  and  some  of  their  characteristics  are  discussed.  In 

3 

Section  -iTT",  the  simulation  program  is  described.  By  using 


the  program,  graphs  and  data  can  be  produced  that  represent 
realizations  of  the  process.  The  program,  which  is  written  in 
BASIC,  is  for  an  HP-85, 


\ 


II.  A  Description  of  the  Process 

The  process,  as  it  is  described  here,  generates  a  random 
track  of  predetermined  length  that  connects  two  predetermined 
end  points.  The  track  originates  at  a  predetermined  starting 
point,  the  first  point,  and  terminates  at  a  predetermined  stopping 
point,  the  final  point.  The  track  is  continuous  between  these 
two  points  and  consists  of  a  sequence  of  straight  line  track 
segments.  Some  representative  tracks  are  shown  in  Figure  1. 

A  sequence  of  rectangular  coordinate  systems  is  used  to 
define  the  process.  The  first  is  a  reference  system  for  the 
first  track  segment,  the  second  is  a  reference  system  for  the 
second  track  segment  and  so  on.  The  x-axis  of  the  first  system 
is  coincident  with  the  line  joining  the  first  and  final  points, 
the  origin  is  at  the  mid-point  of  the  line  and  the  system  is 
oriented  so  that  the  final  point  is  on  its  positive  x-axis. 

For  each  successive  system,  the  x-axis  is  coincident  with  the 
line  joining  the  last  point  of  the  preceding  track  segment  and 
the  final  point,  the  origin  is  at  the  mid-point  of  the  line 
and  the  system  is  oriented  so  that  the  final  point  is  on  its 
positive  x-axis.  Each  track  segment  is  determined  by  two  quan¬ 
tities:  an  angle  0  and  a  length  r.  The  angle  6  is  measured 
clockwise  from  the  direction  of  the  positive  y-axis  of  the 
reference  coordinate  system  of  the  track  segment.  It  is  the 
relative  course  that  will  generate  the  track  segment.  The 
length  r  determines  the  track  segment  length. 
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The  process  is  defined  as  follows:  For  each  track  segment 
except  the  last  a  relative  course  6  is  chosen  from  a  distribu¬ 
tion  with  density  function 

fg (8)  »  (1/2tt)  b/(a  -  c  sin  8) 

where  0  <  8  <  2ir,  2c  is  the  distance  between  the  first  and 

final  points,  2a  is  the  track  length  with  a  >  c  >  0  and 
2  2  k 

b  *  (a  -  c  )  .  Next,  a  length  r  is  chosen  from  a  conditional 
exponential  distribution  with  density  function 

f R | e (r I  9 )  =  [1/6 (8) ]exp[-  r/5 (8)3 

where  0  <  r,  6(8)  =  (a  -  c)5/(a  -  c  sin  8)  and  5  is  the  maximum 
expected  track  segment  length.  Using  the  parameters  ak,  bk  and 
ck  defined  below,  if  r  >  2ak,  then  *  bk  /  ^ ajc”cjc  s*n  9^)  •  Ir* 
this  case,  the  (k+l)st  track  segment  is  the  line  that  joins  the 
last  point  of  the  kth  track  segment  to  the  final  point  and  the 
process  terminates.  If  r  <  2a jc,  then  =  r.  In  this  case,  the 
last  point  of  the  kth  track  segment  and  the  final  point  deter¬ 
mine  the  (k+l)st  reference  coordinate  system  and  the  process 
continues.  The  parameters  in  the  above  expressions  are  defined 
as  follows:  2 ak  is  the  remaining  track  length  with  2a^  =  2a 

and  2a^  »  2ak-i"rk-l  ^or  *  >  ^ck  the  d^-stance  from  the 

last  point  of  the  (k-l)st  track  segment  to  the  final  point  with 

2  2k 

2c^  -  2c  and  b^  ■  (a^  -ck  )  .  The  program  that  is  described  m 
Section  III  can  be  used  to  simulate  this  process. 

Figure  2  illustrates  the  generation  of  a  three  segment 
track.  The  ellipse  on  which  the  second  segment  terminates 


First 

Point 


Figure  2.  A  three  segment  track  and  a  focusing  ellipse 


Figure  3.  The  first  rectangular  coordinate  system,  the 
associated  polar  coordinate  system  and  the 
reference  ellipse. 
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and  the  third  segment  begins  has  one  focus  at  the  last  point  of 

the  first  segment  and  the  other  at  the  final  point.  With  2c2 

the  distance  between  these  two  points,  2a2  *  2a  -  the 

2  2  2 

remaining  track  length  and  b2  »  a2  “  c2  •  t^ie  ellipse  is 

2  2  2  2 

defined  through  x  /a2  +  y  /b2  *  1  in  the  second  reference 

2 

system.  It  is  also  defined  by  r  ■  b2  /(a2  -  c2  sin  0)  in  the 

polar  coordinate  system  with  its  origin  at  the  last  point  of 

the  first  segment  and  its  polar  axis  parallel  to  the  y-axis  of 

the  second  reference  system. 

As  a  basis  for  discussing  the  relative  course  and  track 

length  distributions  that  are  associated  with  the  process, 

consider  the  polar  coordinate  system  with  its  origin  at  the 

first  point  and  its  polar  axis  parallel  to  the  y-axis  of  the 

first  reference  system.  Next,  consider  a  reference  ellipse 

with  one  focus  at  the  first  point  and  the  other  at  the  final 

2 

point  that  is  defined  by  r  =  b  /(a  -  c  sin  9).  The  coordinate 
system  and  the  ellipse  are  illustrated  in  Figure  3. 

The  distribution  of  the  relative  course  0  and  the  condi¬ 
tional  exponential  distribution  of  the  length  R  are  both  deter¬ 
mined  by  o(9)  =  (a-c)/(a  -  c  sin  9),  the  ratio  of  the  value  of 
the  radial  coordinate  of  the  ellipse  at  9  to  its  maximum  value 
at  t/2.  Note  that  when  a  track  segment's  relative  course  is 
tt/2  the  segment  is  directed  toward  the  final  point.  The  dis¬ 
tribution  of  the  relative  course  0  is  determined  by  the  require¬ 
ment  that  the  ratio  of  the  probability  that  a  relative  course 
is  between  9  and  9  +  A9  to  A0  is  proportional  to  p(9).  The 


conditional  exponential  distribution  of  the  length  R  is  deter¬ 
mined  by  the  requirement  that  the  ratio  of  5(8),  the  expected 
track  segment  length  for  8,  to  5,  the  expected  track  segment 
length  for  it/ 2,  equals  p(6). 

If  the  first  point  and  final  point  are  coincident,  then 
c  *  0  and  the  reference  ellipse  becomes  a  circle.  In  this  case, 
the  distribution  of  the  relative  course  0  i  uniform  with  den¬ 
sity  function  fg(9)  *  1/2tt  where  0  <  8  <  2ir  -ne  relative  course 
0  and  the  length  R  are  independent  and  the  ;  ibution  of 
R  is  exponential  with  density  function  fR(r)  "  (1/6) exp  (-r/6) 
where  0  <  r. 

Consider  a  process  that  is  defined  in  the  same  way  except 
that  the  length  of  a  track  segment  is  always  chosen  to  be  equal 
to  the  value  determined  for  R.  For  this  process,  a  is  not  the 
track  length,  there  is  not  a  length  constraint,  and  the  final 
point  is  not  an  end  point,  there  is  not  a  termination  constraint 
The  final  point  is  a  reference  point  of  the  process,  literally, 
focal  point  of  the  process.  If  c  =  0,  the  process  is  identical 
to  the  random  tour  process  described  by  Washburn  in  Reference  1. 
If  c  >  0,  it  is  a  generalization  of  that  process. 

By  inspection  of  the  density  function,  it  can  be  seen  that 
both  the  mean  and  the  mode  of  the  relative  course  3  are  equal 
to  */2.  For  the  constrained  process,  this  implies  that  the 
most  likely  next  to  last  and  last  track  segment  relative  course 
values  are  tt/2  and  -ir/2.  In  Appendix  1,  it  is  shown  that  the 
expected  value  of  the  length  R  is  a6/(a+c) .  For  the  uncon¬ 
strained  process,  this  implies  that  the  average  track  segment 
length  is  a6/(a+c). 


A  new  process  could  be  generated  by  redefining  the  para¬ 
meters  a  and  c.  For  example,  one  could  use  the  definition 
a  *  2c  where  c  is  again  the  distance  between  the  first  and  final 


points,  in  this  case  the  degree  of  focusing  of  the  process 
would  be  constant.  However,  relative  to  the  process  described 
here,  it  is  conjectured  that,  for  track  lengths  greater  than  2a, 
the  process  would  concentrate  track  points  about  the  final  point, 
for  track  lengths  less  than  2a,  the  process  would  concentrate 
track  points  about  the  first  point.  In  the  simulation  program, 
for  a  track  of  m  segments  where  m  >  2,  the  maximum  length  of 
the  first  m-2  track  segments  is  computed.  And,  it  is  also  con 
jectured  that,  for  a  track  length  less  than  2a,  the  probability 
that  the  length  of  the  last  track  segment  would  exceed  this 
statistic  would  be  greater  than  that  for  the  process  described 


! 


III.  The  Simulation  Program 

The  simulation  program  requires  three  inputs:  The  distance 
D  between  the  end  points,  the  track  length  TL  and  the  maximum 
expected  track  segment  length  <5 .  The  program  is  for  an  HP-85 
and  is  written  in  BASIC  for  that  computer.  The  program  listing 
is  in  Appendix  2.  After  loading  the  program  and  initiating  RUN, 
a  user  is  prompted  for  these  inputs  and  advised  of  their  limits. 
An  example  of  the  prompting  with  a  set  of  inputs  is  shown  below. 


DISTANCE  <0105*10)  ? 

19 

TRACK  LENGTH  <0<TL150>‘ 
59 

DELTA 

1 


Table 


1.  The  input  limits  are  based  on  display  and  memory 


requirements 


As  a  simulation  progresses,  the  developing  track  is  displayed 
on  the  HP-85  CRT.  After  a  simulation  is  completed,  a  user  is 
given  the  option  of  displaying  track  data  or  printing  the  track 
and  the  track  data.  In  the  CRT  display  of  a  track,  the  starting 
point  is  on  the  left  side  and  the  stopping  point  is  on  the  right 
side.  In  the  printed  output,  a  track  is  rotated  90°  so  that  the 
starting  point  is  at  the  top.  Figure  4  shows  sections  of  the 
printed  output.  In  the  columned  data,  I  is  the  track  segment 
number,  $  is  the  track  segment  course  for  a  starting  point  that 
is  due  west  of  the  stopping  point,-  R  is  the  track  segment  length, 


i 
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8  is  the  track  segment  relative  course  and  X  and  Y  are  the 
coordinates  of  the  last  point  of  the  track  segment  in  the  first 
coordinate  system.  If  the  starting  point  is  considered  to  be 
due  west  of  the  stopping  point,  then  the  positive  x-axis  of  the 
first  coordinate  system  is  directed  east  and  the  positive  y-axis 
is  directed  north.  Recall  the  origin  of  the  first  coordinate  system 
is  at  the  mid-point  of  the  line  joining  the  starting  and  stopping 
points. 

For  the  track  shown  in  Figure  4,  the  termination  criterion 
was  satisfied  at  the  simulation  of  the  49th  track  segment.  In 
the  first  line  below  the  columned  data,  9  is  the  average  track 
segment  relative  course  where  -90 °<  9  <  270°,  R  is  the  average 
and  RM  is  the  maximum  track  segment  length  for  the  first  49 
track  segments.  In  the  second  line,  4>(51)  is  the  relative  course 
of  the  51st  and  last  track  segment  and  R(51)  is  its  length. 

If  a  track  simulation  reaches  97  segments  without  termina¬ 
tion  and  the  termination  criterion  is  not  satisfied  for  the  98th 
segment,  the  simulation  is  stopped.  In  this  case,  the  values 
on  the  first  line  below  the  columned  data  refer  to  all  of  the 
98  track  segments  and  "TRUNCATED"  is  on  the  second  line.  For 
both  cases,  the  average  relative  course  9,  the  average  track 
length  R  and  the  maximum  track  length  RM  relate  to  the  uncon¬ 
strained  stochastic  process. 

Figures  5,  6  and  7  show  sections  of  the  printed  output  for 
three  additional  simulated  tracks.  The  plots  illustrate  the 
effect  of  the  value  of  6  on  the  plot  scale.  Doubling  the  value 
of  5  doubles  the  distance  between  points  on  the  plot.  This 
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increases  the  likelihood  that  the  HP-85  CRT  will  display  a  com¬ 
plete  track.  The  plots  illustrate,  to  some  extent,  the  increased 
focusing  that  occurs  as  the  value  of  TL  approaches  the  value  of 
D. 

A  second  simulation  program  to  generate  track  data  in  a  non- 
graphical  but  more  efficient  manner  was  written.  A  preliminary 
analysis  of  some  data  generated  by  the  simulation  has  been  made. 
In  Table  2  below,  some  estimates  p  with  associated  95%  confidence 
intervals  are  given  for  p,  the  probability  that  the  last  track 
segment  length  will  exceed  the  maximum  track  segment  length  of 
the  preceding  track  segments. 


D 

TL 

5 

P 

95% 

confidence 

interval 

10 

50 

1 

.60 

.57 

<  P 

<  .63 

10 

50 

2 

.55 

.52 

<  P 

<  .58 

10 

50 

5 

.53 

.50 

<  P 

<  .57 

10 

20 

1 

.44 

.40 

<  P 

<  .47 

10 

20 

2 

.43 

.  39 

<  P 

<  .46 

10 

20 

5 

.45 

.42 

<  P 

<  .48 

Table  2. 

Estimates  p 

for  p. 

the  probability 

that 

.  the 

last 

track  segment  length  will  exceed  the  maximum  track 
segment  length  of  the  preceding  track  segments. 
Each  estimate  is  based  on  a  sample  size  of  900. 
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Figura  4 
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DISTANCE*  5.00 
TRACK  LENGTH*  50.00 
DELTA*  1.00 
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DISTANCE* 10 . 80 
TRACK  LENGTH*  30.00 
DELTA*  2.00 
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Figure  6.  A  program  generated  track. 
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Appendix  1  Some  Mathematical  Relationships 

In  the  simulation,  values  of  9  are  generated  by  using  the 
inverse  transform  method.  This  method  requires  the  use  of 
the  cumulative  distribution  function  FQ (0)  and  it  can  be  found 
as  follows:  For  -tt/2  <  0  <  tt/2  , 

0 

FQ(0)  *  /  (l/2ir)  [b/(a-c  sin  9)]  d9 

0  -tt/2 

*  2[b/ (a+c) ] (tan”^[ (a  tan  (9/2)-c)/b]  -  tan”*[- (a+c) /b] } 

For  it/2  <  0  <  3tt/2,  Fq(9)  -  1  -  Fq(tt-0). 

The  expected  value  of  R  for  the  unmodified  process  that  is 

given  in  Section  II  can  be  found  as  follows: 
tt/2 

E(R)  *  2  /  E(R|e»9)  fQ(0)  d0 

-tt/2 

t/2 

*  /  ( S/t )  b (a-c) [1/ (a-c  sin9)^]  d0 

-V2 

»  (25/tt)  [a/ (a+c)  ]  {tan”*[b/ (a+c)  ]  +  tan'^f (a+c) /b] } 

*  a 5 /a+c 

If  the  termination  criterion  is  satisfied  on  the  first 

track  segment,  the  track  will  consist  of  only  two  track  segments 
2 

With  r  (9)  *  b  /(a-c  sin  0),  the  probability  P,  of  this  event 
can  be  found  as  follows: 


*/2 

P2  -  2  /  P(R  >  re(0) |0  *  0)  f0(9)  de 

-ir/2 

ir/2 

-  2  /  exp[-re(8)/5(9)]  f0(O)  d0 

-tt/2 

»  exp[- (a+c) /<S] 
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Appendix  2.  The  HP~85  Simulation  Program 
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COM  Fa00>,R<te0>,TUd0>,XU 
eo), r<l00> 

OISP  "OISTANCE"  ;  "  <8";CHR*U 
38) i “0" ; CHR#< 163) ; " 10) " ; 
INPUT  08  Cl *0^2 
IF  0> 10  OR  D<0  THEN  28 
OISP  “TRACK  LENGTH"; "  <D";CH 
R*<60>;  ■  TL "  ; CHR# <188)  ;  "S0>"; 
INPUT  L 

IF  L.>50  OR  LOO  THEN  5©  ELSE 
Al* L/Z 

OISP  "DELTA";"  U " ; CHR#< 1 88) 
; CHR#< 18) ; * ) " ; 

INPUT  D1 

IF  O1O0  THEN  38 
RAN00MI2E 

X<0)«-C1  4  Y<0>»0  8  F<0.)«0  8 
R<0)*8  8  T<0>»0 
81*0  8  82*0  8  B3*0  8  Rl*8  8 

A*A1  8  C*C1  8  S*0  8  V*0  8  W* 

0 

GCLEAR  8  E»10*SOR<O1) 

SCALE  -E,E- -< .?3*E>, ,75*E 
MOVE  Cl,  8  8  LABEL  “o" 

MOVE  -Cl,e  8  LABEL  "o" 

FOR  1*1  TO  S3 
U*RND  8  G1*0 
IF  U-s*  3  THEN  220 

U*U-  5  8  Gl *1 

Nl*SQR<<Rl+Cl)/',Al-CU  >  8  N2 
*SQR  •'  1  -C 1  tC  1  •'  <  A 1  *  A 1 ;  > 

T *2*ATN<  N2tTPN<P  I  *U-ATN<.  N 1 )  > 
♦Cl '  At ) 

IF  Gl*l  THEN  T*PI-T 
P«<A1-C1 >'< Al-C1*SIN<T) ) 

02*0 1*P 

R*-^D2tL0G< 1-RNO) > 

F*T  ♦$ 

GOSUB  300 
IF  A2>C2  THEN  340 
R* < A*A-C*C > / 1  A«*C*S  1 N •; T  ,< ) 
G03U8  308 
82*1 

DRAW  X<I),Y<1) 

A*A2  8  C*C2  8  S*-ATN2if , X>  8 
s<  I  ;>*R 

H*T  4  GOSUB  778 
r<i>-H 

H*F  8  GOSUB  778 
F<  I >*H 

IF  62*1  THEN  458 
Rl*HAX<R, R1 > 

V*V*<T-V).'I  4  H*H+<R-ta>/l 
NEXT  I 

OISP  4  OISP  “TRUNCATED"  8  K * 
I  4  83*1  4  GOSUB  330  8  GOTO 

ioa 

ORAN  Cl. -6 

K* 1 ♦ 1  4  R*K>*2*C  8  T  •  fc  > »?0  4 
H»S*PI/2  8  GOSUB  778  4  F  i. K  > 


478  GC'SUB  850 

480  OISP  4  OISP  “OISP  OftTft  ifsH> 

U 

498  INPUT  8* 

590  IF  THEN  520 

510  GOTO  680 

520  OISP  USING  530  ;  CHR#< 15) , CH 
R$( 16> 

530  I WAGE  ,3',X,  "I“,3X,  A,4X,  "P**, 
5X.ft.5X,  "X"  *  6X  .•  MY" 

548  FOR  I»l  TO  K- 1 
550  OISP  USING  560  ;  I..F(I),R(I) 
,  T  <  I  J> ,  X  <  I  >  ,  Y  t  I  > 

560  IMAGE  , / , 20 , 2X • 32 , X , 20  2D , 2X 
, 32, X, 30.20.X, 30. 20 

578  NEXT  I 

580  OISP  USING  590  ,  CHR*0  44;>,V 
,  CHR#t!210>  ,  W,  R1 
590  IMAGE  A, -*“,32,2X,A, 

20.2D.2X, "RH*",20  20 
680  IF  83*1  THEN  648 
610  OISP  USING  620  ;  CHR*«ll5>,K- 
FCK>,K,R<K> 

620  IMAGE  ,  ' ,  A ,  ,  20 ,  *>*•'•  32 , 2X 

,  "R<  20,  “>**•,  20. 20, 5s. 

630  GOTO  668 

640  OISP  USING  658  ;  “  TRUNCATED* 

650  IMAGE  , st  9ft, 9^, 

668  IF  B 1 *8  THEN  680 
670  CRT  IS  1  8  GOTO  ?10 
680  OISP  8  OISP  -PRINT  PLOT  &  OA 
T  ft  • Y- N> " ; 

690  INPUT  84 
700  IF  B**"Y"  THEN  720 
710  ALPHA  8  EMO 
720  PRINT  USING  730  ;  D,L,DJ 
730  IMAGE  . Zs , “DISTANCE*" , DO . DO, 
.'V  ,  "  TRACK  LENGTH*  "  ,  DOG  .  00  •  x/1 
•  "DELTA**,  DO.  00,  3-' 

740  CRT  IS  2  8  61*1 
750  IRAPH  8  COPY 
760  GOTO  520 

770  H-RT0<HJ  8  H*RI1O<'H,360> 

730  IF  h-;0  THEN  H*H+360 
790  RETURN 

300  X  <  I )*X< 1-1  '♦PtSIN < F  > 

310  Y *  I  , «Y(  I- 1  > +R*COS  F  > 

320  X*Cl-X<  I  >  0  V*-Y<fl  ,' 

830  C2*S*JR<  XfcXVfTY  \s2  8  A2*A-R/2 
340  RETURN 

350  H*V  4  H*R  TO  <  H  >  0  U*H 
860  RETURN 
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